Determination of complex absorbing potentials from the electron self-energy 
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The electronic conductance of a molecule making contact to electrodes is determined by the 
couphng of discrete molecular states to the continuum electrode density of states. Interactions 
between bound states and continua can be modeled exactly by using the (energy-dependent) self- 
energy, or approximately by using a complex potential. We discuss the relation between the two 
approaches and give a prescription for using the self-energy to construct an energy- independent, 
non-local, complex potential. We apply our scheme to studying single-electron transmission in an 
atomic chain, obtaining excellent agreement with the exact result. Our approach allows us to treat 
electron-reservoir couplings independent of single electron energies, allowing for the definition of a 
one-body operator suitable for inclusion into correlated electron transport calculations. 

PACS numbers: 03.65.Nk 05.60.Gg 73.63.-b 
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I. INTRODUCTION 

Initiated by experimental advances, interest has been 
growing in the first-principles description of quantum 
transport through nanojunctions formed by a single- 
molecule bridge between electrodes acting as electron 
reservoirs (see Ref. H! for a recent overview). Molecular 
electronic structure, including its response to external 
fields, are well-described by ab initio methods if many- 
electron methods are used. Recent studies also point to 
the need for a detailed quantum chemical approach to 
predict current-voltage characteristics for electron trans- 
port across single molecules^ii^ii^'^. Unfortunately, the 
interaction of the molecular energy states with the bulk 
electrode density of states makes explicit treatment of 
the many-electron problem intractable. 

Similar coupling of bound and continuum states oc- 
curs in many diverse processes in chemistry and physics 
but a completely first principles description is gener- 
ally difficult. For example, with conventional, basis-set- 
dependent methods, it is not feasible to efficiently de- 
scribe both bound and continuum states simultaneously. 
However, it has been known for some time that a reduced 
description for a finite number of degrees of freedom can 
be formally achieved by adding to the uncoupled Hamil- 
tonian of the selected subsystem a non-Hermitian effec- 
tive interactionii^i^ii^. As a result of this extra term, 
the line spectrum of the uncoupled Hamiltonian evolves 
so that the sharp energy levels become resonances with 
broadenings and shifts that depend on the form and 
strength of the coupling to the continuum. 

An exact non-Hermitian interaction can be obtained 



through the Feshbach-FanOi 
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projection operator 



technique and is essentially the self-energy of Green's 
function method o^^'^'^i-'^^'^^i^''''^^ . It is not entirely clear 
how best to proceed in the general case but there are now 
standard numerical methods to extract a first-principles 
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self-energy for bulk electrodes in the single-particle ap- 
proximation^. For certain model systems, analytical ex- 
pressions are availabl o^^i^° . If we could use this coupling 
to the reservoirs, we could describe quantum transport 
through a nanojunction using sophisticated many-body 
methods on a region (e.g., an "extended" molecule in- 
cluding the molecule plus some part of the leads) while 
retaining a single particle description of the bulk elec- 
trodes. This would open up the possibility to improve on 
the commonly employed but controversial single-particle 
description of single-molecule devicea^iii^iS. However, ap- 
plication of one-electron self-energies in many-electron 
calculations is not straightforward, as the single particle 
levels coupled to reservoirs through the self-energy terms 
have no immediate corresponding quantities within cor- 
related electronic theories. Hence we seek a method that 
is able to transmit the information contained within the 
self-energies, e.g., the coupling to electron-reservoirs, di- 
rectly to a many-electron description. 

The formally exact complex-scaling technique pro- 
vides just such an alternative^. By applying a 
complex-coordinate transformation to the Hamiltonian 
via smooth-exterior-scaling the Schrodinger equation is 
transformed; the eigenvalues remain the same but the 
wave functions obey different boundary conditions. This 
yields exact resonant positions and widths. However, for 
the specific applications we have in mind it is opera- 
tionally more convenient to use another general means 
that is commonly employed and only requires the addi- 
tion of an energy- independent, phenomenological com- 
plex potential to the Hamiltonian^i^. 

These potentials are typically local in space and purely 
imaginary, with negative imaginary part; they vanish in- 
side the subsystem region (well-described by the usual 
basis set), and grow rapidly away from that region. The 
negative imaginary walls cause asymptotic damping of 
resonant eigenfunctions, preventing them from extending 
to infinity. The wave function becomes square-integrable, 
obviating the need to describe free-particle states or 
the associated continuum. Hence, standard bound-state 
methods can be applied to open systems. Because these 
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potentials effectively absorb particles that would other- 
wise escape to infinity, they are known as complex ab- 
sorbing potentials (CAPs). 

Calculations may depend rather sensitively on the de- 
tailed form of the CAB2ii2^i2^i2ii^, and one must there- 
fore be careful in constructing and parameterizing the 
complex potential. Semi-classical arguments have been 
used to suggest a form^&iSliSS^ and constraints on the form 
are knownS2i. Parameters are sometimes fit to experi- 
ment; more frequently they are numerically optimized so 
that the stationarity condition for the complex-value of 
the resonant energy is satisfied;^. Deviations occur be- 
cause in practice complex potentials not only absorb par- 
ticles but also cause artificial reflections. In fact though, 
Riss and Meyer— and Moiseye^^ pointed out that there 
is a relation between complex-scaling and CAPs in the 
limit of zero reflections which may be used to introduce 
exact reflection-free CAPs; the parameters in their func- 
tional form are determined by further stability condi- 
tions^. 

While the utility of the complex potential approach has 
been demonstrated many times^i^, we require a technique 
which allows us to avoid a search over parameter space if 
we are to study molecules interacting with large electron 
reservoirs. Further, we would like to make use of the well- 
understood techniques for calculating the first-principles 
self-energy for bulk electrodes. This suggests relating 
the self-energy with an energy-independent complex ab- 
sorbing potential. This CAP could then be used in a 
many-body calculation as an additional one-body poten- 
tial, similar to other applications of complex absorbing 
potentia b--! - ? . This would allow us to apply many-body 
methods on a region while using the self-energy trans- 
formed CAP to couple to the reservoirs. 

In this paper, our purpose is to present a method to 
determine a CAP starting from the self-energy. As as first 
application of our approach we calculate the transmission 
spectrum of independent electrons propagating through 
a model electrode-molecule-electrode junction that has 
been previously reported from calculations using both 
exact methodsi^i^ and a numerically optimized CAP—. 
Remarkable agreement with the exact result is found, 
thus substantiating a link between the exact self-energy 
and an energy-independent CAP. 

The structure of the paper is as follows. In Sec. |TT1 we 
introduce background theoretical material for the cou- 
pling of bound states to a continuum and derive our 
CAP from the self-energy. Section IIIII summarizes the 
application of our method to single-electron transmis- 
sion through a simple model for an electrode-molecule- 
electrode junction described in Sec. IIII Al The results of 
our calculations are presented in Sec. IIII Bl We conclude 
with a summarizing section. 



II. THEORY 

In order to exploit the apparent relation between the 
self-energy and the complex absorbing potential, we as- 
sume that an energy-independent Hamiltonian with the 
"correct" broadened and shifted energy levels and states 
contains all the relevant physics from the interaction to 
the continuum. We must simply define these energy lev- 
els and states and then build the operator. 



A. Coupled States and Broadened Energies from 
the Self-Energy 

Suppose we have a bare Hamiltonian Hq which de- 
scribes a closed, finite system. The fact that the system 
is closed and finite means that Hq will have bound states 
\Xi) with sharp energies e^. That is, we have 

HolXt) e^lXi), (la) 
{Xj\Ho - e,(x,|, (lb) 
iXilXj) = Sij. (Ic) 

If we now couple our system to a continuum (or, equiv- 
alently, give it infinite extent), the energy levels will 
become Ui — ei + Si — iXi , where Si is a shift in the posi- 
tion of the energy level and A; is a width, and the states 
described by \xi) will change. Our goal is to obtain these 
new energy levels and states. 

This can be achieved in a formally exact way using the 
self-energy E from Green's function theory; the proper 
states and energy levels can be obtained by solving the 
Dyson equation. That is, we solved 

[Ho + nUJ^)]\^P^) = (2a) 
[i70 + S(w,)] = ^^{^^\■, (2b) 

the real part of uji gives the position of the i^^ resonance 
including the shift from e^, and the imaginary part gives 
the level broadening. There will almost always be more 
solutions of this equation than there were eigenvalues of 
our original bare Hamiltonian. The extra solutions cor- 
respond to states dominated by the continuum to which 
we have coupled, and in these states we are not partic- 
ularly interested. We need some process, then, to find 
the states and energy levels that best correspond to the 
states and energy levels of the bare Hamiltonian. 
To do this, we solve the related problem 

[i/o + ASK^)] 1^,^) = u^l^^), (3a) 
(0^|[i/o + AS(c.,^)] = u^icj^^l (3b) 

At A = 0, we have the original states \xi) and energy 
levels ei, while at A = 1, we have the target states and 
energy levels. We simply let A go adiabatically from to 
1, thereby obtaining the desired states \4>i) and {(j)i \ and 
energy levels uji. 
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B. Definition of the CAP Operator 

In principle, we would like to build an energy- 
independent complex potential W such that the Hamil- 
tonian Hq + W has the uJi as eigenvalues, and the 
and \4>i) as left- and right-hand eigenvectors. Unfor- 
tunately, we cannot actually construct such a Hamilto- 
nian. This is because if and lipj) are eigenvectors of 
the same Hamiltonian, they must satisfy {(j}i\ipj) — Stj. 
But because we obtain and \ipj) at different ener- 
gies UJi ^ ujj, they do not obey this biorthogonality 
relationship. We must therefore consider constructions 
which yield eigenvalues that are approximately the uJi 
and eigenvectors that are approximately the {(f)i \ and {tpi) 
while obeying the biorthogonality constraint. If we have 
approximate left- and right-eigenvectors ((f>[ \ and and 
eigenvalues w^, we can define 



W = Y,mu;'M\-Ho. 



(4) 



Then by construction, Hq+W will have the desired eigen- 
values and eigenvectors so long as {(p'ili^j) = Sij. We 
consider several ways to proceed. 

The simplest approach is to build Hq + W using the 
eigenvectors of Hq. That is, we could define 



(5) 



This assumes that getting the proper broadening and 
shifts is all that is really needed to capture the proper 
physics. 

Our next step is to introduce the dual spaces to {tpi) 
and {4>i\. That is, we have vectors {Tpi\ and \(f>i) which sat- 
isfy {Tpi\ipj) = Sij and {4>i\(i>j) = Sij. With these in hand, 
we can now build two approximations to W , namely 



i 



(6a) 
(6b) 



Both Ho+W'^ and Hq+W'^ have the correct eigenvalues; 
the former gives the correct right-hand eigenvectors but 
rotated left-hand eigenvectors, and the latter does the 
opposite. Provided we are able to define the dual spaces, 
([Sa|) and (j6b[) give different potentials. However, we ex- 
pect them to yield identical results as the same physics is 
carried by either the left- or right- hand eigenfunctions. 
Finally, we consider a further alternative. We define 



W 



(7) 



which assumes that symmetrizing will have useful effects. 
In this case, however, neither the eigenvalues nor the 
eigenvectors will be correct. 

Note that, to lowest order in GqY., where Gq is the bare 
Green's function (that is, Gq — — HQ + iri]~^), the var- 
ious states we use to build the sundry approximations to 



W are the same, so one expects that all these energy- 
independent approximations to the self-energy would 
yield roughly similar results. Also, we have really defined 
the complex potential in Hilbert space, and not as some 
explicit real-space function. In general, W may be non- 
local and have both real and imaginary parts, like S(i?) 
and unlike the phenomenological complex potentials in 
use. 



C. Matrix Formulation 

It may prove helpful to put everything in matrix lan- 
guage briefly. We suppose we have a bare Hamiltonian 
matrix Hq, which has eigenvectors X so that the eigen- 
value problem of ([1]) becomes 



HoX 
X+Ho 

xtx 



Xe, 

eXt, 

1. 



(8a) 
(8b) 
(8c) 



Once we add the self-energy matrix S(-E), the Dyson 
equation of ^ becomes 



V| [Ho + S(c^,)] = ^.V|. 



(9b) 

We build up total eigenvector matrices U and from 
the individual eigenvectors, and similarly build a total 
eigenvalue matrix u). Note that by taking the transpose 
of (j9bp . we see that if both Hq and S(i^i) are symmetric 
matrices, then V* = U so that V"'' = U"'". 

Our various approximations for the complex potential 
then give us 



Wo 

w 



If Vt = U''^, then W is symmetric. 



Xa;Xt - Ho, 
UwU-i - Ho, 
V-^ujV^ - Ho, 
W'^ -t- W'* 



(10a) 
(10b) 
(10c) 

(lOd) 



III. APPLICATION 

In the previous section, we presented a formal deriva- 
tion of an energy-independent CAP from the self-energy 
that couples the subspace Hamiltonian to the continuum 
states. In what follows, we test the various approxima- 
tions for the complex potential W and examine its struc- 
ture in the case of transmission of electrons through an 
atomic chain. 



A. Model System 

To simplify the calculations and to compare with previ- 
ous results based on the more conventional form of CAPs, 
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we turn to a simple Hiickel model for an atomic chain 
with one orbital per atomic site and nearest neighbour in- 
teractions. Physically, this is intended as a (very simplis- 
tic) treatment of an electrode-molecule-electrode system. 
Typical molecular junctions are made of 7r-conjugate car- 
bon chains bonded by anchor groups to metal electrodes. 
Our model Hamiltonian reads 



H = 



H 



M 



(11) 



where i?L(R) describes the left (right) electrode, Hm 
describes the molecule, and V describes the electrode- 
molecule coupling. Using the ((c^)^) operator that 
creates (annihilates) an electron on the i'^ site of region 
k, the various components of the Hamiltonian are 



R 



oo 

E[£«(cf)tcf-7H(cf)tcf±J, (12b) 



N 



E[eM(cf)tcf -7M(cf)tcf^J, (12c) 



V = -[(ri(cf)tcf + rfl(40tcf) + h.c.].(i2d) 



i=i 



Unless stated otherwise, we take sl = £r = £m = ^o, 
IL = IR = IM = 1, and Tl = Tr = 7/2. The same 
model has been studied using the exact self-energ y ^^i^** , 
and has also been investigated with a complex potential 
of the usual type?'^ , 

The Hamiltonian of has an infinite number of de- 
grees of freedom due to the infinite size of the electrodes. 
We relate to our previous discussion in sections [J and |TT] 
by considering A'l (-/Vr) sites from the left (right) elec- 
trode and the N sites representing the molecule. The 
bare Hamiltonian matrix Hq of this subsystem is thus of 
dimension iVr = A^l + N + N^, and is tridiagonal; the 
diagonal elements are all filled with cq, and the subdiago- 
nal (and superdiagonal) elements are given by —7, except 
for the elements that represent the coupling between the 
molecule and the electrode, which are given by —T. 

The self-energy needed to account for the rest of the 
(infinite) electrodes {i.e., to describe the couphng to the 
continuum) can be calculated exactly for this modeU^. 
The expression reads 



i^ /T^ ) : |77|<1, 
\/rf — 1) : else, 



where 



V 



27 • 



(13) 



(14) 



Note that if we do not explicitly include any of the elec- 
trode sites, i.e., we take Ni^ = or N^i = 0, we must 
scale the self-energy for that electrode by (r/7)^. 



In any event, the total self-energy matrix 1^{E) can 
be written as the sum of the self-energy matrices for the 
left and right electrodes, respectively Tli,{E) and ^^(iJ), 
which are given by 



[^n{E)], 




« = i = 1, 
else, 

i^j^NT, 

else. 



(15a) 



(15b) 



B. Results 



Using the self-energy of pS]) . one can calculate the 
transmission function through the molecule via 



T(i;) -tr(ALGARGt) 



(16) 



where G = [E - (Hq + 'S{E))]-'^ is the Green's function 
matrix, and Al and Ar are the spectral densities for the 
left and right electrode, respectively, defined in terms of 
the self-energies as 



^L(R) 



= i s 



L(R) 



^L(R) 



(17) 



If one uses a complex absorbing potential instead of the 
self-energy, one simply replaces the self-energy matrices 
Sl(£^) and Sr (£') with complex potential matrices Wl 
and Wr when defining the spectral densities (Al and 
Ar) and the Green's function (G). 

We proceed with the analysis of the transmission spec- 
trum for our model. To this end, a useful measure of the 
deviations of the transmission Tw using our approxima- 
tions to W from the exact result Ts is provided via 



AT 



ma.x\Tw{ri) - T^{i])\ 

V 



(18) 



and 



77 is a dimensionless parameter defined in (|14p . 



(19) 



1 . Tests of Different Approximations to the Complex 
Potential 

We have several possibilities to consider. The first 
thing we should point out is that, while with the usual 
local complex potentials it is trivial to define the poten- 
tial in each electrode, this is not quite as obvious with 
our potentials. That is, if we follow our procedure as 
outlined in Sec. |TT1 building the total complex potential 
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W directly from Hq + Sl + , we find that it is not, 
in general, local, and we cannot uniquely define Wl and 
Wr. For the moment, if we use our procedure to build 
Ho + W in one step, we write 











Wl = 





















( 







Wr = 


















On the other hand, we could simply use Hq + 5]l(r) 
to build Ho + Wl(r'), in which case we have a clear 
and unambiguous way of building the potential due to 
each electrode at the cost of doing twice as much work. 
We wish to consider all four possibilities for building W 
given T:{E) as listed in ([TU)) . and for each we wish to see 
whether it suffices to build W directly and from it obtain 
Wl and Wr via ([20| or if we must build Wl and Wr 
separately from the beginning. 

In Fig. [TJ we show results from each of our four approx- 
imations (see Eg. [T0|). wherein we generate W directly; 
Fig. [2] is identical except that we build Wl(r) separately 
through Sl(r). In all cases, we use TVl = A'r = 100 and 
iV = 12, and compare to the exact result generated by 
using the self-energy. 

Let us first examine the results obtained by construct- 
ing W directly and extracting Wl and Wr according to 
(1^0)) . Although none of these results are perfect, we note 
that most of them are fairly reasonable. Using the un- 
perturbed eigenvectors to build W is clearly inadequate 
to describe the transmission, but already at this crud- 
est level of approximation, we see that we can obtain 
some qualitative features. Unsurprisingly, simply using 
the correct eigenvalues is sufficient to put the resonance 
peaks in the right positions and produce a fair amount of 
broadening. The unperturbed eigenvectors cannot, how- 
ever, account for a quantitative description at energies 
far from the resonance peaks, as they yield overestimated 
broadenings. 

In order to do a reasonable job over the whole en- 
ergy range, we apparently need to use the correct states. 
Using the exact states to build W'^ and W*^ gives us 
a transmission function that, except for some relatively 
small oscillations which are most apparent near the min- 
ima, is really quite good. W''' and W*^ yield identical 
results, as anticipated, and excluding some effects at the 
band edges that we discuss separately, our deviation mea- 
sures read AT^^ax = 0.058 and ATave = 0.014. Av- 
eraging the two, to yield a symmetric complex poten- 
tial performs almost as well with ATmax = 0.071 and 
ATaye = 0.025. But Hq -I- W has neither the correct 
eigenvalues nor the correct eigenstates, again pointing to 
the importance of using the right states. 

Though we can already get reasonable transmission 
functions by building W directly so long as we use the 
proper states and energy levels, there are still noticeable 
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FIG. 1: Transmission functions with A^l = A^r ~ 100 and 
= 12 from various approximations used to build W di- 
rectly. The exact result is shown with the heavy line, while 
the result from our complex potential is given in the dashed 
line. We recall that rj — {E — eo)/2'y. 



errors which we would like to eliminate. More important, 
however, is that we have no a priori way of defining the 
potential due to each electrode separately, and our use of 
(PO)) to separate Wl from Wr is largely arbitrary. We 
can solve both of these problems by building Wl and 
Wr separately and combining them to build W at the 
end. 

Note first that we must exercise some care if we follow 
this course. If we were to simply use the eigenvectors 
X of the unperturbed Hamiltonian to do this, we would 
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FIG. 2: Transmission functions with A'l = Nb. ~ 100 and 
= 12 from various approximations used to build Wl and 
Wr, separately. The exact result is shown with the heavy 
line, while the result from our complex potential is given in 
the dashed line (for W") or in dots (for W''' , W", and W). 
We recaU that rj = {E ~ eo)/27. 



obtain Wl = Wr — 1/2 W, which is clearly nonsensical 
and is why the transmission function built from W*' is 
meaningless. To amplify on this point, note that under 
this construction, the spectral densities Al and Ar are 
diagonal in the same basis as is the Green's function. If 
the eigenvalues of Hq + Wl and Hq + Wr are cui = 



ei + 1/2 5i — 1/2 iXi, then in the diagonal basis we have^l 

(At,'),-,: — A,; 



{Gk^ = 



1 



E ~ El - Si ~ iXi ' 
1 



E ~ e, ~ S, + iX, 
The transmission function is thus 

A? 



[E - e, - + A? 



(21a) 
(21b) 

(21c) 
(21d) 

(22) 



In other words, it is the sum of A^t independent 
Lorentzian resonances; this explains the jagged nature of 
the calculated transmission function, as well as the fact 
that T{E) is far too large. This should be contrasted 
with the exact result and with our results using other 
constructions of W, in which interference between states 
suppresses all but N peaks in the transmission function. 

When, however, we use W' or W^, the results are 
for all practical purposes perfect; ATmax — 0.006 and 



0.001. While we see a slight degradation in 



quality from using W instead (with ATmax — 0.093 and 
ATave = 0.009), any of these three approaches would 
be reasonable to take. Once again, the importance of 
generating the right states is clear. 

In the remainder of this work, we will use only W"^, 
built by constructing Wl and Wr separately, as the re- 
sults clearly indicate that this is the best way to construct 
a complex potential among all those that we have con- 
sidered as demonstrated in the middle panel of Fig.[2P^. 



2. Effects of Electrode Size 

One of the most interesting features of our complex 
potential is that the results for the transmission func- 
tion show an unexpected insensitivity to TVl and TYr . In 
Fig. [31 we show results for iV = 12 and for varying elec- 
trode sizes TVl and TVr. The results are virtually indistin- 
guishable; we can even choose to use TVl = A^r = and 
still obtain a remarkably accurate transmission function. 
Of course, the exact transmission function is strictly inde- 
pendent of these parameters, but this observation points 
to the redundancy of an absorption grid outside the area 
of interest in similar applications. With a complex po- 
tential of the more usual form, this is not necessarily the 
case. In common implementations, a finite absorption 
region, whose extent depends on the precise form of the 
CAP, is essential for the damping of the resonant wave 
function. 

To some extent, this success is illusory. Even our op- 
timal complex potentials do not work as well far from 
resonance or close to the band edges as they do near a 
resonance; any residual deviations from the exact results 
occur at these energies. While the small errors observed 
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Electrode Atoms 
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FIG. 3: Transmission as a function of rj = {E — eo) for TV — 
solid line; our results are the dots. 

may not be of practical relevance, the electrode extent 
still plays the major role if stringent convergence is re- 
quired. We demonstrate this with the following exam- 
ples. 

Resonance Minima - A careful perusal of Fig. [3] reveals 
that our results are not quite as good near the minima 
of the transmission function as they are in the resonance 
peaks. This effect can be seen more readily by decreas- 
ing N, which decreases the number of resonances and 
increases the range of energies far from any of the reso- 
nances. In Fig. 21 we show results for A'^ = 3 and varying 
TVl and Ar . We can now see that while once again we can 
use A^L = Nyi = to describe the transmission function 
at resonance, we must use rather large leads if we are to 
obtain accurate results far from any of the peaks. Indeed, 
even with A'l = A^r = 60, there are still some residual 
oscillations about the exact result. On the other hand, 
the deviation measures ([18]) and p9)) decrease mono- 
tonically from AT^^ax = 0.037 and ATave = 0.019 to 
ATmax = 0.016 and ATave — 0.004 as the electrode size 
varies from A^l = A^r = to A^l = A^r, = 60. 

Band Edges - One way to see the difficulty in describ- 
ing the transmission function near the band edge is to 
consider the eigenvalues of the matrix T{E) given by 

T(i;) = AlGArG^, (23) 

whose trace yields the transmission function. The ex- 
act matrix derived from the self-energy has exactly one 
non-zero eigenvalue at each energy; this eigenvalue is of 
course numerically equal to the transmission function at 
that energy. This is to be compared with the results from 
our calculations using W''' , displayed in Fig. [51 Near the 
band edges, we have two non-zero eigenvalues at each en- 
ergy, and their sum gives us the calculated transmission 
function. As the energy moves away from the band edge, 
the smaller of the two eigenvalues goes to zero and we 
approach the exact result. The precise energy range that 
marks convergence within a given tolerance depends on 
the size of the included electrodes. This is readily seen 
by the eigenvalue plots in Fig. [51 as A'l and A'r increase. 
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12 with varying TVl and A'^r. The exact result is given by the 



Ideal Wire - A stricter test of our complex potential can 
be done by calculating the transmission function when 
the electrode-molecule coupling parameter F is set equal 
to the other coupling parameter, 7. In this case, the 
exact result is that T{E) = 1, as it should be for an 
ideal wire with no potential sources for scattering. We 
expect that our complex potential would find this case 
more difficult to describe, and that the sensitivity to A^l 
and A'r (or, as we cannot distinguish molecule sites from 
electrode sites, the sensitivity to A^t) should be much 
larger than in our previously considered cases. We show 
results for A^t = 72 and for A^t = 212 in Fig.H Note the 
reduced range of the ordinate in these plots. In either 
case, we see small oscillations around the exact result, 
and we note that the results are considerably worse near 
the band edge than they are near the middle of the band. 
The results improve by increasing A^x, again pointing 
to the role of the electrode size. We believe that the 
problems near the band edge might be related to the van 
Hove divergence in the electrode density of states at i? = 
±27. It is worth noting that these same difficulties in 
describing transmission through an ideal wire are found 
when phenomenological complex potentials are used in 
the calculation—. 



3. Structure of the Complex Potential: Non-locality 

Now that we have seen how well our complex potential 
works, it is worthwhile investigating its structure. We 
will continue to focus on W''' , and on Wl and Wr built 
separately. 

With the bare Hamiltonian and the self-energy that 
we are using, Wl and Wr take a particularly simple 
structure. For Wl, only the first row is non-zero, and 
for Wr only the last row is non-zero; for symmetry rea- 
sons it is of course the case that Wr can be obtained 
from Wl by left-right reflection (that is, reflection about 
the skew diagonal) . Elements of the first row of Wl for 
TVl = A^R = 30 and A^ = 12 are shown in Fig [3 Note 
that they show a repeating pattern in which each imagi- 
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30 Electrode Atoms 60 Electrode Atoms 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 



FIG. 4: Transmission as a function of ?; = (i5 — eo)/27 for A'' = 3 with varying A^l and A'r. The exact result is given by the 
solid line; our results are the dots. 

Electrode Atoms 20 Electrode Atoms 




FIG. 5: Transmission function and non-zero eigenvalues of T(77) for N — 12 and varying A^l and Nr. The exact transmission 
function (Ts) is displayed in the heavy line; the transmission function calculated using the complex potential (Tw) is displayed 
in the light line, and the two non-zero eigenvalues of T(j7) are given by the dotted line. We recall that rj — {E ~ eo)/2"f. 



nary element is followed by a real element, and each real 
element is followed by an imaginary element. Further, 
while the magnitude of the element typically decreases 
as one moves from left to right, this is not always the 
case and there are some relatively significant elements 
coupling one end of the left electrode to the other end 
of the right electrode. Our complex potential, in other 
words, is strongly non-local. 

The structure of our complex potential is in sharp dis- 
tinction to the more common complex potentials. For 



example, the potential used by Kopf and Saalfrank^S, has 
the usual structure, namely, the only non-zero elements 
are negative imaginary and on the diagonal. They van- 
ish within the molecule and the opposite electrode, and 
increase in magnitude as one moves away from the molec- 
ular region. 

We have seen that our best complex potentials are non- 
local, but that the elements far from the diagonal are 
rather small. This naturally raises the question of just 
how important the non-locality is. We can study this 
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FIG. 6: Transmission as a function oi rj = {E 
T{ri) = 1. Notice the range of the ordinate. 



eo)/27 with F = 7 for two different values of A't. The exact resuh is that 



Matrix Elements of First Row of Wi 




FIG. 7: Elements of the first row of Wl for Nl ^ Nr = 30 
and A'^ = 12. Some elements are real; the others are imagi- 
nary. None have both real and imaginary parts. 



issue by simply constructing the full, non-local complex 
potential, but, once it is built, removing all but the di- 
agonal (i.e., the fully local part) and some number Ns 
of the subdiagonals. As Ns approaches iVr — 1, then, we 
approach the fully non-local CAP in a smooth way. In 
Fig. [SI we display results for the transmission function 
with TVl = A^R = 30, A^ = 12, and several values of Ng. 

Clearly, the local approximation is wholly inadequate, 
and while including only a few subdiagonals improves 
the results dramatically, there are still deviations from 
the exact result even with Ns — 30. Since there are no 
such deviations for the fully non-local complex potential, 
it is clear that the strong non- locality plays a critical role 
in mimicking the effects of the self-energy. 



IV. DISCUSSION AND CONCLUSIONS 

We have given a prescription for calculating a com- 
plex absorbing potential from the self-energy, thus en- 
abling one to describe coupling to continuum states in 
an energy-independent way. In our test system, the non- 
empirical complex potential gives reasonable results even 



at the simplest level of approximation, but by obtaining 
the potential for each electrode individually and using the 
proper states to construct the potential, we can build a 
non-local complex potential that yields for most practical 
purposes essentially the exact transmission function. 

We have several main observations to make from this 
initial work. 

First of all, it is essential that the complex potential 
yield not only the correct energy levels but also the cor- 
rect states. This is not surprising, as one function of the 
complex potential is to force free-particle states to look 
rather like bound states outside of the area of interest. It 
is no coincidence, then, that inside the area of interest, 
the states should be greatly modified from the original 
bound states to closely resemble the resonant eigenfunc- 
tions. 

Secondly, in cases such as the one we have examined, 
where there are multiple self-energies in the problem, we 
should construct a different complex potential for each 
self-energy, and combine them to form a total complex 
potential at the end. While we can get away with a 
somewhat ad hoc construction if we have some physical 
insight, the results are nevertheless not as good as when 
we build the complex potential for each self-energy sep- 
arately. 

Third, we can, surprisingly, eliminate most of the ab- 
sorption grid outside the subsystem region of primary 
interest. This was shown with the explicit elimination of 
all electrode atoms in our simple system, and it is appar- 
ently another consequence of using the correct resonance 
structure and states. The main effect of increasing the 
number of explicit electrode atoms (which presumably 
decreases the burden on the complex potential) is to im- 
prove the description at energies off resonance or near 
the band edge. The great utility of such a significant re- 
duction in the number of required degrees of freedom is 
readily appreciated if one considers that the bulk metal 
electrodes can be removed almost completely from elec- 
tron transport calculations with CAPs. 

Our complex potential works only imperfectly for the 
ideal wire. Again, most of the problems are at the band 



10 



Subdiagonals 



1 Subdiagonal 





2 Subdiagonals 



5 Subdiagonals 



1 

0.8 
0.6 
0.4 
0.2 





-0.5 

n 



0.5 1 



1 

0.8 
0.6 
0.4 
0.2 





-0.5 
11 



0.5 1 



10 Subdiagonals 



30 Subdiagonals 





FIG. 8: Transmission function with ^Vl — A^r = 30 and A'^ = 12 for various values of Ns- The exact result is shown with the 
heavy hne, while the result from our complex potential in a semi-local approximation is given in the lighter line. We recall that 
V = {E~ eo)/27. 



edge, but if we do not include a sufficient number of sites, 
the description even in the middle of the band has some 
unphysical oscillations. Of course, in such a system, there 
are no resonances and one would not really expect a com- 
plex potential to describe the physics perfectly anyway. 
Nevertheless, convergence of the results can be tested as 
one increases the number of sites explicitly treated. 

Finally, our complex potential is strongly non-local, 
and that non-locality is essential in its functioning. A 
semi-local approximation is not terrible, and works quite 
well on resonance, but the full non-locality is needed if 
we are to accurately describe energies off resonance. 



All of this seems to suggest that, so long as one is care- 
ful in how one uses a self-energy transformed complex 
potential of the sort we have introduced, one should eas- 
ily be able to obtain results within the required accuracy 
without having to introduce parameters. 

A limitation to our approach is the requisite prior 
knowledge of a self-energy at some level of theory. There- 
fore, our method is ill-suited to applications for which 
the self-energy is difficult to construct. However, there 
is a wide range of problems that can be considered. In 
the simple case that many-body effects are of major con- 
cern only within a subsystem'*'^ while the surrounding 
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environment can be adequately treated at the single- 
particle level, our self-energy transformed complex po- 
tential should be easy to construct and would allow us to 
readily describe the embedded subsystem at a high level 
of theory. In particular, for quantum transport across 
molecular junctions it has been shown that many-particle 
effects^ can play a critical role, but including the contin- 
uum of the reservoirs has proven challenging. Our poten- 
tials should enable us to treat the extended molecule with 
an accurate wave function approach while still including 



a reasonable description of the bulk electrodes. 



Acknowledgment 

We would like to thank Science Foundation Ireland 
(SFI) for funding this work. One of us (E.H.) was funded 
through the SFI UREKA program. 



^ G. Cuniberti, G. Fagas, and K. Richter, eds., Introducing 
Molecular Electronics, vol. 680 (Springer, Berlin, 2005). 

2 R Delaney and J. C. Greer, Phys. Rev. Lett. 93, 036805 
(2004). 

^ G. Fagas, P. Delaney, and J. C. Greer, Phys. Rev. B 73, 
241314 (2006). 

* M. Koentopp, K. Burke, and F. Evers, Phys. Rev. B 73, 
121403 (2006). 

^ M. Albrecht, B. Song, and A. Schnurpfeil, J. Appl. Phys. 

100, 013702 (2006). 
" K. S. Thygesen and A. Rubio, cond-mat/0609223. 

R. Santra and L. S. Cederbaum, Phys. Rep. 368, 1 (2002). 

* J. G. Muga, J. P. Palao, B. Navarro, and I. L. Egusquiza, 
Phys. Rep. 395, 357 (2004). 

^ N. Moiseyev, Phys. Rep. 302, 212 (1998). 
^° W. Domcke, Phys. Rep. 208, 97 (1991). 
" H. Feshbach, Ann. Phys. (N.Y.) 5, 357 (1958). 
12 H. Feshbach, Ann. Phys. (N.Y.) 19, 287 (1962). 
1^ U. Fano, Phys. Rev. 124, 1866 (1961). 

J. S. Bell and E. J. Squires, Phys. Rev. Lett. 3, 96 (1959). 

L. Fonda and R. G. Newton, Ann. Phys. (N.Y.) 10, 490 

(1960). 

1^ L. S. Cederbaum, Phys. Rev. Lett. 85, 3072 (2000). 
1^ R. Santra and L. S. Cederbaum, J. Chem. Phys. 117, 5511 
(2002). 

1* S. Feuerbacher, T. Sommerfeld, R. Santra, and L. S. Ceder- 
baum, J. Chem. Phys. 118, 6188 (2003). 

1^ G. Cuniberti, G. Fagas, and K. Richter, Chem. Phys. 281, 
465 (2002). 

2° V. Mujica and M. A. Ratner, Chem. Phys. 326, 197 (2006). 

21 T. Seideman and W. H. Miller, J. Chem. Phys. 96, 4412 
(1992). 

22 D. Neuhauser, J. Chem. Phys. 103, 8513 (1995). 

2^ J. G. Muga and B. Navarro, Chem. Phys. Lett. 390, 454 
(2004). 



2^ D. Macias, S. Brouard, and J. G. Muga, Chem. Phys. Lett. 

228, 672 (1994). 
25 U. V. Riss and H.-D. Meyer, J. Phys. B 28, 1475 (1995). 
2** A. Vibok and G. G. Balint-Kurti, J. Chem. Phys. 96, 7615 

(1992). 

2^ D. E. Manolopoulos, J. Chem. Phys. 117, 9552 (2002). 
2* B. Poirier and J. T. Carririgton, J. Chem. Phys. 118, 17 

(2003) . 

2" U. V. Riss and H.-D. Meyer, J. Phys. B 26, 4503 (1993). 
^° U. V. Riss and H.-D. Meyer, J. Phys. B 31, 2279 (1998). 
^1 N. Moiseyev, J. Phys. B 31, 1431 (1998). 
^2 R. Santra, L. S. Cederbaum, and H. D. Meyer, Chem. 
Phys. Lett. 303, 413 (1999). 

Y. Sajeev, M. Sindelka, and N. Moiseyev, Chem. Phys. 
(2006). 

G. Fagas, G. Cuniberti, and K. Richter, Phys. Rev. B 63, 
045416 (2001). 

^5 A. Kopf and P. Saalfrank, Chem. Phys. Lett. 386, 17 

(2004) . 

'^^ If I '5) is an eigenvector of a Hermitian operator, then its 
Hermitian adjoint 1^')^ will be the same as its dual 
Because Ho + E is non-Hermitian, the dual of one of its 
eigenvectors will not be the same as its adjoint. In an abuse 
of notation, we write — even if l^*) comes from a 
non-Hermitian operator. That is, the bra is always meant 
to indicate the Hermitian adjoint of the ket, rather than 
its dual. 

'^'^ That there are no factors of one half in the Green's func- 
tion is because each of Wl and Wr contributes 1/2 5^ — 
1/2 i A,. 

We could of course use instead, as it is equally good. 
A. Kopf, private communication. 

T. M. Henderson, J. Chem. Phys. 125, 014105 (2006). 



